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8 : Abstract 
(N 

^ ' Recent experimental data on the radiative decays B V'j, where y is a light vector meson, 

■ find small isospin violation in i? ^ K*'y while isospin asymmetries in B ^ fry are of order 20%, 

. with large uncertainties. Using Soft-Collinear Effective Theory, we calculate isospin asymmetries in 

these radiative B decays up to ©(l/m^), also including 0{vas) contributions from nonperturbative 
charming penguins (NPCP). In the absence of NPCP contributions, the theoretical predictions for 
Q^' the asymmetries are a few percent or less. Including the NPCP can significantly increase the isospin 

^. asymmetries for both B Vj modes. We also consider the effect of the NPCP on the branching 

O I ratio and CP asymmetries in B^ V^'y. 
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The rare radiative B decays, B V'y, where V is a hght vector meson, are important in 
heavy flavor physics because the dominant processes are due to the flavor changing neutral 
current. Isospin asymmetries are interesting observables for testing the Standard Model 
(SM) and investigating new physics in the flavor sector. The isospin asymmetries for B — »• 
K*^ and B ^ p'j are deflned to be 

K*%) - T{B- K*--f) 
T(B° + T{B- K*--f) ' 

2T(b' PO7) + T{B- ^ ■ 

In these asymmetries, the decay rates are averaged over charge conjugate modes. Recent 
experimental measurements flnd [l| 

A^' = 0.03 ± 0.04, A(;_ = 0.26 ± 0.14, (2) 

where the average values for the decay rates for B K*{p)'y are taken from the Heavy 
Flavor Averaging Group The isospin asymmetry for B i^*7 is consistent with zero 
within an error of a few percent. The data suggests that the asymmetry in i? ^ p7 
is signiflcantly larger, but because of large uncertainties it is not yet possible to draw a 
deflnitive conclusion. The work in this paper is motivated by the question of whether an 
anomalously large isospin asymmetry in i? — p7 can be understood within the SM. In 
particular, we calculate subleading contributions to the leading QCD factorization theorems 
for B —>■ V'-f to see if they can explain the observed asymmetries. 

In the heavy quark limit, the leading amplitudes are factorizable However, 
isospin violating asymmetries come from 0(1/ rrib) suppressed power corrections for which the 
factorization is necessarily more complicated. In this paper we calculate 0{l/mh) corrections 
to the asymmetries using Soft-Collinear Effective Theory (SCET) ^], which provides a 
systematic power counting. In addition, possible endpoint divergences in these higher order 
corrections can be regulated without imposing an arbitrary infrared cutoff by including the 
zero-bin subtraction of Ref . . Previous QCD analyses of isospin asymmetries in radiative B 
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FIG. 1: Nonperturbative charming penguin (NPCP) contributions for (a) B — > Vj and (b) B 
M1M2 arise when x = 1 — x « 4m^/m^, in which case the charm quark pair is in the threshold 
region. 



2 



decays appear in Refs. 10|, lUl, 112] • The main difference between our analysis and previous 
work is the inclusion of nonperturbative charming penguin (NPCP) contributions, which 



JJ, UM ■ (For 



are already known to play an important role in nonleptonic B decays \lc 
alternative point of view, see Ref. fl^). 

The NPCP contributions to B ^ V'j are depicted in Fig. [T]-(a). In certain kinematic 
regimes, the invariant mass of the charm quark pair in the loop in Fig. [I]-(a) is near the 
threshold, 2mc, in which case the charm quark pair is described by nonrelativistic QCD 
(NRQCD) and additional interactions need to be taken into account. As pointed out in 
Ref. [13], contributions from this regime are suppressed by only vas{2mc) compared to the 
leading contribution. Here v is the relative velocity of the charm quarks in the threshold 
region. Therefore, the NPCP contribution to the isospin asymmetry could dominate over 
other I /nib suppressed contributions. In this paper, we calculate the isospin asymmetries 
including the NPCP along with l/rrih suppressed contributions. We also calculate the NPCP 
contributions to the branching ratio and CP asymmetries for B^ V^j. 

In the absence of NPCP contributions, the theoretical predictions for A^* and Aq_ are 
no larger than a few percent. Including the NPCP contributions can significantly modify 
the theoretical predictions for the isospin asymmetries. We will see below that the NPCP 
contribution can be factorized using SCET, and the result expressed in terms of nonper- 
turbative matrix elements. The NPCP matrix elements are fitted to available data on the 
isospin asymmetries, A^* and Aq_, the CP asymmetry for B^ p^7 (recently measured 
by Belle 17|), and the branching ratio for B^ [2]. The predictions for A^* and Aq_ 

are of order 10%, with uncertainties large enough that both predictions are consistent with 
experiment. However, the NPCP does not predict a large difference between A^ and Aq_, 
as suggested by the central values in Eq. ([2j). 

The isospin asymmetry in — > V'y can arise from either the mass difference of the 
spectator quark in the B meson or the electric charge difference when the spectator quark 
emits the photon in the final state. However the isospin asymmetry due to the mass difference 
is negligible because it is (9((m„ — rriii) / Aqcb) and therefore of order 1% or smaller. So the 
dominant piece comes from electromagnetic (EM) interactions with the spectator quark. 

In order to describe the isospin breaking corrections to B ^ V'y, we need the following 
effective weak Hamiltonian 



w 



G p 
71 







,p=u,c 



(3) 



where the operators are 



Oip = {pb)v^A{qp)v-A 
03,5 



O 



77 



iqb)v~A 

q'=u,d,s,c. 



I'q')vtA, 



ga^^F^,(l + 75)&, 



04,6 

(9) 



{Pfsba)v-A 

grrib 



V-A, 



E ( 

q'=u,d,s,c,b 



Here q is the ci or s quark, the CKM factor is A^^^ = VpbV* s^YidV ±A = 7'^(1 ± 75). 



(4) 
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FIG. 2: Various isospin breaking contributions in full QCD. Here crosses denote another possible 
photon emissions from the spectator quark. If we do not consider the long distance contribution 
in the diagram (c), all the contributions are power-suppressed by 0(A/mft). 



For the asymmetric contributions, the photon radiates from either the initial or final 
spectator quark as shown in Fig. [2l If the photon is radiated from the initial spectator 
quark, we need SCET operators of the type O^^''^'^^ = Xn^by^rJ^^n where Cn(Xn) is n{n)- 
coUinear field^ and n and n are light-cone vectors satisfying = = 0, n ■ n = 2. 



The analysis of factorization for these operators appears in Refs. 19|, |20] and the Wilson 
coefficients at next-to-leading order (NLO) in as have been calculated in Refs. 18|, |20 . 



The leading operator in SCET only contributes to longitudinally polarized vector meson 
production, but in i? — V'y the vector meson must be transversely polarized. Transversely 
polarized vector mesons can be produced from subleading operators that are higher order 
in the SCET expansion parameter, A. The relevant effective weak Hamiltonian in SCET is 



(1,49) 
W,SCET 



^ p i=l ^ 



(5) 



where O^^'^"^^ are 



(l,4g) 



o. 



{l,4g) 

2,3 

(l,4g) 



(6) 



q—u,d,s 



o 



(l,4g) 
5,6 



q'=u,d,s 

Here the superscript '1' denotes suppression by one power of A compared to the leading 
operator, Wn{n) is a collinear Wilson line in n(n)-direction, and 1^ is a ultrasoft (usoft) 
Wilson line. The subscript outside the square brackets denotes that a delta function which 
fixes the momentum fraction, x, is included in the bilinear operator: 



^^^Wjix- 



V 



2E,, 



(7) 



Our conventions are the same as Ref. 



18[ except that we have exchanged n ^ n compared to that paper. 
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where V = n ■ V is a derivative operator taking the largest momentum component and Ey 
is the energy of the produced vector meson. Here, C,n is the power-suppressed component in 
the spin projection of = (^^/4)^„ + {j^'i^/A)^^, where g„ is the colhnear quark field. Using 
the equation of the motion, can be expressed in terms of ^„ as 

i^=Wn^Wllp^^i^. (8) 

The Wilson coefficients in Eq. ([5]) are the same as the leading Wilson coefficients 



Cf of Ref. [ISj] due to reparameterization invariance [21[. Because the isospin-breaking 
contributions from -ffp^/'scET suppressed by Aqcd/w^?) compared to the leading decay 
amplitude, we will suppress EM penguins with C7,8,9,io at tree level and keep only Ci,2,8g at 
one loop order in for our phenomenological analysis. 

The four-quark operators in the effective weak SCET Hamiltonian, Eq. (I5l), contribute 
to the radiative weak decay when the photon is emitted from the initial spectator quark, 
as shown in Fig. [3]- (a). To calculate their contribution, we need to take the time-ordered 
products of o'f''^''^ with the electromagnetic interaction term ^^j^, 

4m = e,g,,y^4(xW4x| + h.c, (9) 

where is a photon field and is the electric charge of the quark. The time-ordered 
products are performed in SCETj with possible offshellness ~ m^A and then matched 
onto SCETii, which describes dynamics with fiuctuations of ~ A^. In the matching, 
rl-collinear fields having large offshellness of O(mfeA) must be integrated out, giving a jet 
function 

(0|T{w4x^(^),XwW^^(0)}|0) = ih{z^)5\z^) j ^e-'=-+/V„.,^(fc_), (10) 

where = n ■ k, k_ = n ■ k, and the momentum of the photon, p'^ = n ■ p^W /2 = nihn^ /2. 
At the lowest order in a^, the jet function is simply Jn-p^ = 

The n- and n-coUinear degrees of freedom are decoupled at leading order in SCET, and 
by using a field redefinition it is possible to decouple usoft degrees of freedom from both. 
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Therefore the n-coUinear piece of the matrix element describing the production of the hght 
vector meson is decoupled from the rT-coUinear and usoft parts. Thus, we can compute B 
to 7 from -ff^i'^scET' which only depends on usoft and n-collinear physics, independently of 
the light-meson production process, which depends on n-collinear physics. After a brief 
calculation, we find 



d'z{^iel)\T{^^W^^i{l-^,)YXiO),C^ilmB), 

i^fBlTlBieT + ^^±(^lu) j dl-Jn-p-X-^-)'pBi^-), 



(11) 



where e'^ = e^'^^"nMryl2 setting e'^i^a _ _]^_ 0^ is a light-cone distribution amplitude 
(LCDA) of B meson 22]. We use the convention of Ref. 23| with n and n interchanged. 

The light meson production is described by the n-collinear part of the matrix elements. 
The matrix elements for production of transversely polarized vector mesons in SCET are 



(12) 



d 



{VAV1)\ i^Wnl'^^.Wlin + irWnl'^l.Wli^ |0) = - ^^^^5^ ^1.^^ (^) , (13) 



where (7^'"^ are chiral-even, twist-3 LCDAs 24] whose asymptotic forms are g'^f'{x) = 3(x^ + 
x^)/2 and g^^\x) = 6xx, where x = 1 — x. 

Combining Eqs. ffTTl) . f|T2l) . and f|T3l) . the matrix element of O^'^'^'' for P 7? for 

example, is 



X 



g^l\x,n) 



9 (a) 



Adx 



9± {x,^m) 



(14) 



where the renormalization scales are roughly given by /i ~ mftAqcD and (Xq ~ Aqcd, and 
Al = e*j_- — is^x^ is the polarization factor for the left-handed p~ and 7. In the case 
of S — i> p7, it is necessary to include one-loop corrections to the jet function, see Ref. [2^ 
for details. 

When the photon radiates from a final state quark (i.e., from the crosses in Fig. [2]), the 
intermediate quark line is hard with offshellness of order m^. Matching onto SCET we 
obtain localized four-quark operators with photons, shown in Fig. [3]-(b). In this case, the 
effective weak Hamiltonian that contributes to the decay amplitudes is 



(1,4<?7) 
W^.SCET 



V ^ n i=l Jo 



1(1.497), 



X,/i), 



where the five-particle operators O^^^'^'^"'^ are 



^ eg/g'r„{^,^}(l + 75 

q'=u,d,s 



1 



(15) 



(16) 
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In Eq. (fTSjl . the Wilson coefficients at NLO are 



where G{sp, x) is 



C5 J 2Ci 



1 + ln- 



1 + ln 



- -G{sp,x) 



-G{sp,x) 



2C. 



' xniB 



(17) 



G{sp,x) = — 4 / dzzz\n{sp 
Jo 



zzx 



ie), 



(18) 



p _ nip/m'^. Similar to Bf, we neglect EM penguins at tree level and only keep Ci_2,8g 



and s 

at one loo p. T he sum of the terms in Eq. f|T7|) proportional to Ci differs by a factor of 3/4 
from Ref . lOl . 

Again, the B —>■ 'j piece of the matrix element factors from the light-meson production 
piece. The n-collinear part in Eq. (fT6l) . describing the vector meson production, gives a 
leading twist LCDA, 4>±{x), whose asymptotic form is 6xx. It can be obtained from the 
following projection: 



LWjix- 



As an example, the B 



2Ev 



J b 



/ twist~2 



-JyEviil 



2Jba N 



■<P±{x). (19) 



K*'^ matrix elements from q'^^'^'^'^^ are 



-*o 



-*o 



fBfk*^BAL{el,r]l) 



(f)±{x,fi) 



(20) 



2 - - - X 

Next we turn to the contributions from NPCP, shown in Fig. [T]-(a). The size of this 
contribution is 0{vas{2mc)) [l5| and therefore is suppressed only logarithmically in the large 
TTLc limit, compared to the power suppression of previously considered 0{l/mfj) contributions. 
Numerically, Aqcd/'^6 and vas{2mc) are roughly the same size, so a priori it is sensible 
to include them at the same order. In fact we will see below that the NPCP gives the 
dominant contribution to isospin violation in B ^ ^7. When x is close to 4s^, long-distance 
interactions govern the charm quark pair in the loop and hence it cannot be separated 
from the B meson. However, the n-collinear piece of can still be decoupled, with the 
dominant part obtained from the leading twist projection of Eq. (fT9l) . The factorization 
process is similar to the treatment in Ref. ISj and we refer the reader to that paper for 
details. The NPCP contribution to the decay amplitude is 



^A('^)(\47|CiOie|i?)NPCP 



G, 



^i'^fv^BVl, / dx6{x - 4sl)M^)Hc-c{x,mB){7m\B), 



ANV2 

where Hcs = 4:CiTcas/{xm%) at lowest order and the six-quark operator O, 
nonrelativistic charm quark fields, c±v, is defined as 

d'y c.^Y^^lT-Y^c^M ^(wMl^ll^Yil - 75)TX^c 



CO 



(21) 



including 



^ cc 



(22) 
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In order to integrate out n-collinear fields with fluctuations greater than Aqcd, we consider 
time-ordered products of O^c with C^l^^, 

t-c = ^V1,| d'z{^{el)\T{Om,4LmB) (23) 



dz^dk 
Air 

where we employed Eq. f fTOj) to obtain the second line of Eq. fl23l) and Occq'b is 

0,,,>,iel,vl,z^,y) = 7„^y^(^,y^,^)|/l7;[^l|7^(l-75)T'^y„c_.(0) (24) 

xc_vF„7.^T"F4c+.(y)c+v7,(l - ^MO). 

The matrix element of Occq'b in Eq. fl23|) is purely nonperturbative. It can be decomposed 
into left- and right-handed polarized contributions, 

j A(0|acv;>(el,r7l,^+,y)|5)e^'^-^+/2 ^ f^^3^ J dl_e-'i-^+/^ (25) 

X [A^(el,r^l)F,^-(fc_,/_) +A«(el,r/l)Fj(fc_,Z_)], 

where Aji = e*j_- + is^^^ is the polarization factor for decay into right-handed flnal 
states. An important point is that the NPCP can give a right-handed polarized contribution 
which is not 0{l/mi,) suppressed. As pointed out in Ref. [26|, any other right-handed 
polarized contributions to the decay amplitude should be suppressed by l/m;,. Therefore 
NPCP could give the dominant contribution to the right-handed polarized decay amplitudes. 
Since the right-handed polarized contribution from NPCP cannot interfere with the leading 
order amplitude which produces only left-handed flnal states, the right-handed contribution 
does not enter into the asymmetries until higher orders. Therefore, we neglect any possible 
right-handed contribution from NPCP in our calculations of the asymmetries. 
Combining Eqs.([2l]), ([23D, and ([25]), we obtain 

= -^\^^^e,,fBf^mlA,{el, r^l)^ dx'^5{x - Asl)Hc,{x), (26) 

where q = d ot s, q' is the B meson spectator quark, and Hcc = xm\Hcc/ {Airas) = Ci H — • . 
Here we have deflned A^^ to be 

A,-,^ = -/ rf/_^^e-('=-+'-)^-/V..,,(fc_)F,i(fc_,/_) (27) 

= - j rf/_J,.,^(-/_)F,^,(-/_,/_)~ j dL^^^tLlIzl, 

Following Ref. [15|], we can power count the size of this correction. The NPCP contri- 
bution is suppressed relative to the leading order term by vas{2mc), and thus of order 
rubV a s i2m c) / Aqcd relative to the other isospin breaking terms considered. Based on this 
power counting, we expect that 

mB0j_(4s^) vrub 



Acc 4s2 Aqcd ' 



(28) 
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The factor 0±(4s^) / (4s^) is formally 0(1) in the power counting of Ref. 15|], but numerically 
0_l(4s^)/(4s^) ~ 4.3, so we keep this factor in estimating m^/Acc. Taking v mi,/ Aqcb ~ 3 we 
find niB/Acc ~ 0.7. The extracted values of mB/Acs are consistent with this naive estimate 
but smaller. For these values of m^/Acg, the NPCP gives significant contributions to the 
isospin asymmetries. 

Finally, there is another interesting isospin-breaking source, a double photon contribution 
with the EM penguin Oy^. It is only available for the decay with an unflavored vector meson, 
i.e., B — > p^'-f. The largest contributions are depicted in Fig. HI Concentrating first on Fig. H]- 
(a), the off-shell photon coming from O^^ produces the vector meson and then an additional 
photon is emitted from the B meson spectator quark. Integrating out the hard photon, we 
can match onto the SCETj operator 

Cr.Or, — / dxC,,ix) ^^W^^l^l + 1,)yX (29) 

where C^^ is equal to Cj^ at tree level. Next we integrate out the n-coUinear fields in the 
time-ordered product with and match onto SCETn. Applying Eqs. ffTOl) . ffTTl) . and 
(PI), we find 

(^±7|C7,07,|i?)2,(a) = ^/B/y^^07,A4el,r^l) / d/_ J„.,,(-/_)0+(/_), (30) 
and, in case oi B p'^'j, 

{pUCr,Or,\B'h,^.^ = ^^(^" ' '^M,^C,,A,{el, vl) [ rf/. J..,,(-/-)0^(/_), 

V2 Zir my J 

(31) 

where Qg = Cg/e with Qu = 2/3 and Qd = —1/3. 

At the lowest order in as, the contribution of Fig. |4]-(b) is the same as Fig. Il]-(a) with n 
and n exchanged. So the result can be written as 

/ 1^ n llf \ QdiQu - Qd) ea mym^ *\ f ji t ( i \ 

{pl_l\C7l07^\B )2T,(b) = 7= -TTlBlp Cj^Aiit^^T]^) / dl+Jn.p^{-l+)(j)B{l+)- 

V 2 ZTT my J 

(32) 

The double-photon contribution is suppressed by a, but enhanced by a factor of ml /my 
due to the virtual photon. Compared to the other isospin-breaking breaking contributions, 
such as those in Eqs. (IT^ and (|2Ui) . this contribution is rather small. 
The isospin asymmetries in Eq. ([1]) are given by 

.V Re{b^-b^) + RRe{b^-b:) 

where 

&r = -^, &r = 9-, (34) 

CyLiv J^y 
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FIG. 4: Leading double photon contribution to the isospin asymmetry, where <8) represents 0^^. 
They are suppressed by am^/my compared to other usual isospin breaking contributions. 



are the leading isospin breaking corrections to the decay amplitude, Ly are the leading 
isospin symmetric decay amplitudes, and cy = 1 for K* , cy = — 1/v^ for p. In Eq. fl33l) . 
b^^ are the corresponding ratios for the charge conjugate modes, and R = |Lyp/|Lyp. L 
can be written as 



V 



e 



(35) 



Lp = 



p 



p=u,c 



where the transition form factor for B V, Cj^, is defined as 



± 5 



(36) 



and we will use = 0.36 ± 0.07 and = 0.27 ± 0.05 [ll| for the numerical analysis. The 
Wilson coefficients y in Eq. (!35l) are 



47r 

Cf fpfyrnB 
N ml 



(37) 



dL 



^ 



+ Y^(a;, Sp) + 



where the hard functions Af\ Gi^si ^"^^ H are available in Refs. , and we followed 

the conventions of Ref. 
Finally we obtain 



h^' = O 

^1^..c /-K 



B 



K - Q '^'^^ 



nib 



XBrrib 



2q 



2 /i _^ fpUhLKP 
■ mb Xsmb 



(38) 
(39) 



Here A^^ = / dl(t)-%{l)/l and we use the following model for </>^(/) [28j 



Apl 



2(^B-1), I] 
m 



irXeiP + p^)\-P + p^ 



(40) 
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where the parameters Xb and ctb are Xb = 460 ±110 MeV, as = 1-4 ± 0.4 at /i = 1 GeV. 
Ki 2 can be written as 



K* 



^2q 



dx—_ 
^ 



x) r 1 
2 



A1{x) +A2{x) 



nas mB 



S{x-4sl)] 



-1 



x) r 1 



p=u,c 



^ 
1 



^2. = / 
Jo 



p=u,c 



(v) 

91 



d (a) 

91 



Adx 



N A„_ 

(x) -Ab / rf/_0+(/_)Jp^(-/_) 



a m{,mB 
TT ml 



(41) 
(42) 



-SpcC.^'^Six-Asl)}, (43) 



(44) 



Here we include only the tree-level contributions to Cabibbo-suppressed terms with in 
Eq. (H2|) because it is numerically comparable to the other term with B^. 

In the convolutions of Ai and (j)±{x)/x in Eqs. ( HTpSil and B4 and g^f^ in Eqs. 



there are endpoint divergences, which can be eliminated with the zero-bin subtractions [9| 




X 

dx 



91' {x) 



X 



^^ 0±(X)+X01(1) 

xP' 

X 



(45) 



The zero-bin subtraction removes infrared divergences from the x integrals. We have dropped 
all finite terms including logarithms associated with rapidity scale dependence. We estimate 
the uncertainty associated with this procedure to be 50%. 

For numerical estimates of b^*''' in Eqs. ( !38l) and ( l39l) . we use the following set of 
parameters: {mb,mc,mB,mK*,mp, fs, fx*, fp, ff* , f±} = {4.8, 1.3, 5.28, 0.894, 0.775, 0.2 ± 
0.03,0.218,0.209,0.175 ± 0.025,0.150 ± 0.025} GeV. The CKM parameters are p = 
0.221±0.064 and fj = 0.340±0.045. All Wilson coefficients and hard functions are evaluated 
at the scale fi = m^, we do not include any renormalization group evolution, and we use the 
asymptotic forms for the vector meson wave function, 
isospin asymmetries in the absence of NPCP contributions are 



and g^l^'^"'\ Our estimates for the 



A 



K* 
0- 



0.04 ±0.02, 



^0- 



0.02 ±0.02, 



(46) 



where the dominant errors come from Xb, C±, and CKM factors. These estimates for A^_ 



are comparable to previous theoretical results [10|, lUJ, l29|]. Comparing to Eq. ([2]), we see 
that A^* is consistent, but Ag_ disagrees by about 1.7 a. 

Next we include the NPCP contribution in our calculation. In addition to the isospin 



asymmetries, the NPCP can contribute to the CP violating asymmetry, A^J__ ITj], and to 
the branching ratio, Br[i?+ p+7] In order to obtain values of uib/J^cc that are not 
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Exp. 


w/o NPCP 


w/ NPCP 




0.03 ± 0.04 


0.04 ±0.02 


0.10 ±0.05 




0.26 ±0.14 


0.02 ±0.02 


0.10 ±0.06 


K- 


0.11 ±0.33 


0.08 ±0.02 


0.07 ±0.13 


Bi[B+ p+7] X 10^ 


0.96 ± 0.24 


1.80 ±0.69 


1.63 ±0.67 



TABLE I: Theoretical predictions with and without NPCP compared to experimental data. 



inconsistent with measurements of these quantities, we perform a least squares fit to all four 
observables. In our calculations of A^J__ and Br [5+ — > p^j], we include only the leading 
order and NPCP contributions, without any 0{l/mb) corrections. An analysis that includes 
the 0(1/ nil,) corrections to all four observables is clearly required but is beyond the scope 
of this paper. 

The results of the fits along with experimental results are shown in Table [B The first 
column lists the observables considered and the second gives their measured values including 
errors. In the third column, we show the theoretical prediction in the absence of the NPCP 
contribution for the values of the parameters given earlier. The last column gives the results 
of the fit with the NPCP included. We extract 



Re 



Arr 



-0.102 ±0.063, Im 



A. 



cc 



0.022 ± 0.255 . (47) 



The of the predictions in column three of Table [T] is 15.2. Including the NPCP, the 
is 12.1, so the overall agreement between experiment and theory is slightly improved. 
Note that after including the NPCP the theoretical prediction for Aq_ increases so that the 
la error band of the experimental result and the theoretical result now overlap. However, 
the prediction for A^* is now significantly increased. The trend suggested by the central 
values of the experimental data, namely a large value of Ag_ and small value of A^*, does 
not seem to be naturally accommodated by including NPCP contribution. However, once 
theoretical and experimental uncertainties are taken into account, the theoretical predictions 
are consistent with both isospin asymmetries. For the range of values of me/Acg obtained 
in our fit, the NPCP does not have significant impact on the theoretical predictions for 
A^|__ and Br [5+ P^l]- Finally, inclusion of NPCP contributions substantially increases 
the uncertainty in all theoretical predictions because the parameter ms/Acg is not well 
constrained. 

These results indicate that the NPCP can increase the isospin violating asymmetry Aq_ 
to bring theoretical predictions closer to current data, while maintaining consistency with 
the other observed asymmetries. It is not possible to obtain predictions for the two isospin 
asymmetries that are in agreement with the central values of both A^* and Aq_. How- 
ever, the uncertainties in the current measurements of all asymmetries are large and better 
measurements are need to determine whether the NPCP is an important contribution to 
B ^ p'-f isospin and CP violating asymmetries. 
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To summarize, we have used SCET to calculate the isospin asymmetries in B — V'y 
decays, including all 0(1/ mh) contributions as well as the 0{vas) NPCP contribution. As 
in nonleptonic B decays |30[], our analysis allows for large NPCP contributions, which could 
account for the large isospin asymmetries measured (with large errors) in B ^ pj. If the 
isospin asymmetries are large and the NPCP is the source of these asymmetries, then we 
also expect the NPCP to contribute to CP asymmetries in B ^ p'-f and give larger than 
expected contributions to the right-handed polarized decay rates in B ^ V'j. We speculate 
that the NPCP could also be responsible for the recently measured enhanced transversely 



polarized decay amplitude for B — > VV 31|- More precise experiments will be needed to 



determine the exact size of NPCP and confirm these predictions. 
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